
CLAIMS 

1. A digital simulation method of a non-linear interaction 
between an excitation source and a wave in a resonator, by 
means of digital signal calculation tools based on equations the 

5 solution of which corresponds to the physical event of a 
phenomenon to be simulated which can translated, at each time 
and at each point of the resonator, by a linear relation between 
two variables representative of the effect and of the cause of 
said phenomenon to be simulated, a method wherein the 

10 impedance or admittance equation is transcribed directly into a 
digital model form enabling to realise a non-linear interaction 
between the two variables of the impedance or admittance 
relation. 

2. A method according to claim 1, for simulating an 
15 oscillating phenomenon, wherein the digital model comprises, on 

the one hand at least one linear part (3), representing the input 
impedance or admittance of the resonator, and, on the other 
hand, one non-linear part (2) modelization the role of the 
excitation source (22) of the phenomenon to be simulated. 

20 3. A simulation method according to claims 1 or 2, for real 

time digital synthesis of an oscillating phenomenon, .wherein, 
from a system of equations between at least two variables 
representative of the behaviour of a wave in the resonator, an 
expression of the input impedance or admittance of the resonator 

25 is established in the form of a linear filter including delays, 
without any decomposition into two-way waves, in order to 
realise at least one linear part (3) of the model. 

4. A method according to claim 3, wherein the linear part 
(3) of the model is coupled with one non-linear part (2) involving 

30 the evolution of the non-linearity as expressed between the two 
variables of the input impedance or admittance relation of the 
resonator. 

5. A method according to claim 4, wherein the linear part 
(3) of the digital simulation model of the impedance or 

35 admittance equation in based on to two elementary waveguides 




fulfilling a transfer function between the two variables of the 
impedance or admittance relation. 

6. A method according to claim 5, wherein the linear part 
(3) with two waveguides of the model is coupled with a loop 

5 connecting the output to the input of said linear part (3) and 
comprising a function (21) involving the non-linearity as 
expressed physically. 

7. A method according to claim 6, wherein the model is 
driven by at least two parameters representative of the non- 

10 linear physical interaction between the source and the resonator, 
by means of a loop connecting the output to the input of the 
linear part (3) and comprising a non-linear function (21) playing 
the part of an excitation source for the resonator. 

8. A method according to one of the previous claims, for 
is real-time sound synthesis of a musical instrument comprising, at 

least, one excitation source with non-linear characteristics and a 
linear resonator, the sound produced by the instrument resulting 
from a coupling, between the excitation source and the 
resonator, expressed at least by a linear impedance or 

20 admittance relation and a non-linear relation between two 
physical variables representative of the effect and of the cause 
of the sound produced, a method where the sound produced by 
the instrument is simulated, in real time, by modelization the 
physical phenomena governing the operation of the instrument, 

25 wherein, to realise said physical modelization, the impedance or 
admittance linear relation is expressed directly and digitally 
between two physical variables representative of the cause and 
of the effect of the phenomenon to be simulated and said 
impedance or admittance relation in digital form is associated 

30 with the non-linear relation between the same variables. 

9. A method according to claim 8 for synthesis of the sound 
of an instrument with complex resonator, wherein the resonator 
is decomposed into a series of successive elements and wherein 
the impedance or admittance relations corresponding 

35 respectively to each element of the resonator are calculated and 




combined in order to obtain a global impedance corresponding to 
the geometry of the resonator. 

10. A method according to claim 8 or claim 9, wherein, for 
real-time synthesis of the sound produced by a wind instrument, 

5 the two variables of the impedance relation are the acoustic 
pressure (p e ) and flow (u e ) at the input of the resonator. 

11. A method according to claim 10, wherein, for an open- 
ended cylindrical resonator, the linear part (3) of the digital 
transcription model of the impedance equation is the sum of two 

10 elementary waveguides having as excitation source the flow (u e ) 
at the input of the resonator, and fulfils the transfer function : 

7 (co ) = Pe(co) = 1 exp(-2ik(co)L) 

eK } Ue(co) 1 + exp(-2ik(co)L) 1 + exp(-2ik(o)L) 

15 wherein: 

- co is the angular frequency of the wave ; 

- Ze(co) is the input impedance of the resonator, 

- Pe(oo) and Ue(co) are the Fourier transforms of the 
dimensionless values of the pressure and of the flow at 

20 the input of the resonator ; 

- k(co) is a function of the angular frequency which 
depends on the phenomenon to be simulated ; 

- L is the length of the resonator. 

12. A method according to claim 11, wherein each of the 
25 two waveguides involves a filter having as a transfer function: 

- F(co) 2 = - exp(-2ik(co)L) 
and representing a two-way travel of a wave, with a sign 
change at the open end of the resonator, each waveguide 
30 corresponding to a term of the impedance equation. 

13. A method according to claim 12, wherein the model is 
driven by the length (L) of the resonator and at least two 
parameters (C, y) representative of the non-linear physical 
interaction between the pressure (p e ) and the flow (u e ) at the 
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input of the resonator, by means of a loop connecting the output 
to the input of the linear part (3) and comprising a non-linear 
function (21) as an excitation source for the resonator. 

14. A method according to claim 11, wherein the non-linear 
function has the pressure and the displacement of the vibration 
formation member as input parameters, and is controlled by at 
least two parameters simulating a player playing. 

15. A method according to claim 14, wherein the playing 
parameters for controlling the non-linear function are: 

- a parameter q characteristic of the mouthpiece and of 
the action of the player on the vibration formation 
member, 

- a parameter y representative of the pressure applied to 
the vibration formation member. 

16. A method according to one of the claims 10 to 15, 
wherein, for real-time synthesis of the sound to be simulated, a 
formulation is realised in the time domain of the angular 
frequency response of the impedance of the resonator, by 
approximation of the losses represented by the filter by means of 
an approximated digital filter. 

17. A method according to claim 16, wherein, to express 
the angular frequency response of the impedance of the 
resonator, a one pole digital filter is used of the form: 



wherein : 

CO 

w = - 

e , fe being the sampling frequency, 

c is the pure delay corresponding to an away or 
return travel of the wave in the resonator, 
- the coefficients bo and a1 are expressed in relation to the 
physical parameters so that F(con = F(tn) for a value 



F(w) = 



b Q exp(-2ixnD) 
1 - exp(-im) 



(13) 



col of the 
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angular frequency corresponding to the fundamental playing 
frequency and another value co2 corresponding to a harmonic, 
and the following differential equation is derived: 

p e (n) = u e (n) - a,u e (n - 1) - b 0 u e (n - 2D) + a.p^n - 1) - b 0 p e (n - 2D) (16) 

18. A method according to claim 17, wherein the 
coefficients bo and a1 are obtained by solving the equation 
system: 

IF^) 2 ! 2 ^ + a 2 - 2a 1 cos(w,)) = b 2 



|F(co 2 ) 2 | (1 + a 2 - 2a 1 cos(tu 2 )) = b 



2 



with |F(co) | = exp(-2acj— L) , said coefficients being given by 
the formulae: 



a, = 



A1 - A2 - V(A, - A 2 ) 2 - (F, ^FJ 
R -F 2 



b 0 = 



_ V2R,F 2 ( Cl - c 2 )(A, - A 2 - /(A, - A 2 ) 2 - (F, - F 2 ) 2 ) 

F,-F 2 



wherein 

c, = cos(tn 1 ),c 2 = cos(tn 2 ),F, = \F(ay,) 2 \ 2 ,F 2 = |F(q> 2 ) 2 | 2 , 
A 1 — F^c^, A 2 = F 2 C 2 

19. A method according to claim 18, for simulating a 
cylindrical resonator instrument, from a physical modelization 
governed by the system of equations : 



15 
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1 d 2 x(t) q r dx(t) /4X , /4X 
'+-^ — ^ + x(t) = ±p e (t) 



w 2 dt 2 (o r dt 
(with the sign + for a reed and the sign - for the lips) 



( I :3/2 \ 

P e (co) - i tan — - - — acco 1/2 L U e (co) 

v c 2 ; 

5 

u e (t) = 1(1 - sign(y - x(t) - 1))sign( Y - p e (t)K(1 - y + x(t))J|y - p e (t)| 

wherein cor is the resonance frequency and qr is the quality 
factor of the reed or of the lips, wherein that said system of 
equations is solved in the time domain from an equivalent 
io sampled formulation of the angular frequency response of the 
displacement of the reed or of the lips and of the impedance 
relation which is translated by the system of equations : 



x(n)=biap e (n-1 )+a 1a x(n-1 )+a 2a x(n-2) (1 8) 

Pe(n)=u e (n)-a 1 u 0 (n-1)-b 0 Ue(n-2D)+aip e (n-1)boPe(n-2D) (19) 

u e (n) = 1 (1 - sign( Y - x(n) - 1))sign( Y - p.(n)Kn)- Y + x(n))^Y - P.(n)| (20) 

said equations being used sequentially by grouping the 
20 terms not depending on the time sample n, in order to calculate 
in succession : 

x(n)=biap e (n-1 )+ai a x(n-1 )+a 2a x(n-2) (21 ) 

25 V=-ai u e (n-1 )-b 0 u e (n-2D)+ai p e (n-1 )-b 0 p e (n-2D) (22) 

W = 1(1 - sign( Y - x(n) - 1))Q) - Y + x(n» (23) 

u e (n) = 1 sign(Y - VX-bCoW 2 + W^bCoW) 2 + 4|y - V| (24) 
Pe(n)=b 0 c 0 Ue(n)+V (25) 

30 



54 



20. A method according to claim 19, for more realistic 
simulation of the produced sound, wherein, by neglecting the 
radiation, the external pressure is expressed as the time 
derivation of the outgoing flow, in the form: 

5 

P-(t) = ^(p.(t) + u.(t)) (26) 

and is calculated, at each sampled time (n), by differential 
between the sums of the internal pressure pe and of the flow ue, 
respectively at the time (n) and at time (n-1). 
10 21. A method according to claim 19 for simulating a 

multimode reed instrument, wherein the calculation of the 
acoustic pressure and of the flow at the mouthpiece is performed 
by sequential resolution of a system of equations wherein the 
displacement of the reed at each time(n) is in the form: 

15 

x(n)=b a1 p e (n-1)+b a2 Pe(n-2)+b aD1 pe(n-D a -1) +a a1 x(n-1) +a a2 x(n-2)+a aD x(n-D a )+a aD1 x(n-D a -1) 
the coefficients aa1, aa2, aaD2, aaD1 being defined by: 

'e *e 'e 

and the coefficients ba1, ba2, baD1 by : 



. C - Ca a . - Cb a 

b a i= — , b 3 2= -, b a D1- a 



fe' fe fe 



l(o r q r et C= J± V 3 
2 (VA7 - 1KA, 



25 P=-co r q r et C= r * (28) 



2 1 

by setting : A n = , A 2 =-« r q r and A 3 =A 2 A 1 q r , 

G> r Vq r 2 +4 2 

the following equations being the same as for a single 
mode reed. 

22. A method according to one of the claims 19 to 21, 
30 wherein that a model is prepared for a cylindrical resonator with 



terminal impedance from the basic model corresponding to a 
cylindrical resonator and being the sum of two waveguides 
involving each a filter having as a transfer function -F(co) 2 =- 
exp(-2ik(co)L), while replacing the expression exp(-2ik(co)L) by 
5 the expression 



Zc being the characteristic impedance 

-^VetZ s the output impedance 

7i R U s (co) 

23. A method according to one of the claims 19 to 22, 
10 wherein, from the impedance model for cylindrical resonator 

instrument and the associated differential equations, other more 
complex impedance models are built for simulating oscillating 
phenomena produced by a resonator of any shape by combining 
impedance elements in parallel or in series and by using digital 
is approximations for an explicit use of the physical variables 
involved in the production of said oscillating phenomena and a 
more flexible control of the result of the simulation. 

24. A method according to one of the claims 19 to 23, 
wherein that, from the basic model of a cylindrical resonator 

20 wherein the angular frequency response of the displacement of 
the reed or of the lips is translated by a system of differential 
equations providing, at each time(n), the displacement x(n) the 
pressure pe(n) and the flow ue(n) at the input of the resonator, a 
model for a conical resonator is built wherein the equation of the 

25 pressure is in the form: 

p e (n) = bCoUe(n) +bciU e (n-1)+bC2Ue(n-2)+bc D Ue(n-2D)+bCDiU e (n-2D-1) 

+ acip e (n-1)+ac 2 Pe(n-2)+acDPe(n-2D)+acDiPe(n-2D-1) (33) 



R(co)exp(-2ik(co)L), wherein R(co) = 



Z c -Z s (co) 



wherein the coefficients bcO, bd, bc2, bcD and bcD1 are 



30 defined by : 



bc 0 = — , bci= - 
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and the coefficients aci, ac 2l ac D and ac D i are defined by : 



_ «.o,+o. ^ 



b ° G ™ ac -h 



10 



by noting : G p = 1 + 



1 



2f e ^ 
e c 



and G m = 1- 



1 



2f -T 



25. A method according to one of the claims 19 to 23, 
wherein that, from the basic model for a cylindrical resonator, a 
model for a short resonator having a length I is built, by an 
approximation of the impedance according to the expression: 



Z,(co)= i tan (k(co)l) = G(co) + icoH(co) 



(34) 



1 - exp 



wherein G(co) = 



" ac if ' 



1 + exp 



- ac if 



and H (co)=-(1 - G(co)). 



26. A method according to claims 23, 24, 25 in 
is combination, for simulating a wind instrument, wherein the 
mouthpiece or the bill is modelled by a Helmholtz resonator 
comprising a hemispheric cavity coupled with a short cylindrical 
pipe and a main resonator with a conical pipe, the input 
impedance of the resonator assembly which may be expressed 
20 as: 

_L_ 



Z e (co) = 



. V 

ICO + 

pc 2 



1 



ICO- 



ICO 



1 + 



itan(k 2 (co)L 2 ) 
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wherein V=— is the volume of the hemispheric cavity, L1 
6 

is the length of the short pipe, L2 is the length of the conical 
pipe, Z1 and Z2 are the characteristic impedances of both pipes 
which depend on their radii, k1(co) and k2(co) take into account 
5 the losses and the radius R1 and R2 of each pipe, and that, from 
the basic model for cylindrical resonator, from its extensions to 
the conical pipe and to the short pipe, a resonator model is 
prepared by expressing the pressure at the mouthpiece or in the 
bill by the differential equation : 



27. A method according to one of the claims 1 to 9, for the 
simulation of an oscillating phenomenon wherein both physical 

15 variables of the linear relation are the strength applied to one 
point of a mechanical system such as a string generating 
vibrations and the speed at this point, wherein that the 
admittance is expressed, at this point, in the form of a 
combination of the admittances of each part of the string, on 

20 both sides of said point, each mechanical admittance being 
obtained from the basic model describing the acoustic 
impedance of a cylindrical pipe resonator, by expressing the 
speed at the point considered of the string in relation to the 
strength applied to this point, where the filter F(co) of the basic 

25 model may be expressed from a bending wave propagation 
model in a having a stiffeness. 

28. A digital device to implement the method according to 
one of the claims 8 to 27, for the simulation of a musical 
instrument generating a sound resulting from a coupling between 

30 a linear resonator and an excitation source with a non-linear 
characteristic, which can be expressed at least by a linear 
impedance or admittance relation and a non-linear relation 
between two variables representative of the effect and of the 



10 



Pe(n)=£bc k u e (n - k) + £bc Dk u e (n - k - 2D) 



(36) 




+ Z ac kPe(n - k) + X> c DkPe(ri - k - 2D) 




cause of the produced sound, said simulation device comprising 
a control element I including at least one gestural sensor 1 
transforming the actions of a player into control parameters, a 
modelization element II including one non-linear part (2) 

5 associated with a linear part (3) and an element creating the 
synthesised sound III, wherein that the linear part (3) includes a 
computing block (31) driven by the length (L) of the resonator, 
having as an input parameter a signal representative of one of 
the variables, cause or effect, computed by the non-linear part 

10 (2), and the transfer function of which is the input impedance or 
admittance of the resonator, in that the non-linear part (2) 
implements a non-linear function (21) driven by at least two 
control parameters and having as input parameters a signal 
representative of the other variable, cause or effect, computed 

15 by the linear part (3) and a signal modelization the role of the 
excitation source, the linear part (3) being thus coupled with the 
non-linear part (2) in a closed loop, and in that the element 
creating the sound III computes a sound signal from signals 
representative of the cause and of the effect of the sound to be 

20 simulated, emitted respectively by the linear part (3) and the 
non-linear part (2) 



